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The results related to non-equilibrium phenomena at the very early and late stages of the 
processes of A+A collisions are presented. A good description of the hadron momentum 
spectra as well as pion and kaon interferometry data at RHIC is reached within the realistic 
dynamical picture of A+A collisions: HydroKinetic Model (HKM). The model accumulates 
the following features: not too early thermalization time; r > 1 fm/c; a developing of the pre- 
thermal transverse flows; the effectively more hard, than in the case of chemical equilibrium, 
equation of state of expanding chemically non-equilibrated multi-hadronic gas; a continuous 
non-equilibrated emission of hadrons. All these factors lead to a good description of the 
mentioned RHIC data, in particular, the observed R ut/Rside ratios, solving, therefore, the 
HBT puzzle in detailed realistic model. 
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I. INTRODUCTION 



The typical spacial and temporal scales of 
the interaction processes in proton-proton colli- 
sions are less or of order 1 fm and 1 fm/c corre- 
spondingly. This is the result of the femtoscopy 
analysis and it is agreed with a simple theoreti- 
cal estimates. The similar experimental measure- 
ments reveals a typical life-time of the system cre- 
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ated in A+A collisions to be at least one order of 
the magnitude larger. The femtoscopy, or inten- 
sity interferometry measurements are associated 
with a set of the points of particle last scattering. 
Therefore the above results mean that, even if one 
tries to consider the collision of nuclei as a some 
kind of superposition of the individual collisions 
of nucleons of nuclei, the secondary hadrons, pro- 
duced in these "elementary" collisions, continue 
to interact during the time interval which is much 
larger than the time-scale related to the individ- 
ual nucleon-nucleon collisions. Then one can con- 
clude that in ultrarelativistic A+A collisions we 
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have to face a phenomenon of the space-time evo- 
lution of the strongly interacting matter produced 
in such processes, and can rise the question as for 
a nature of this matter at the early collision stage. 
The other crucial points are: whether this matter 
becomes thermal in the processes of A+A colli- 
sions, and if yes, how does it evolves and also how 
are observed particle momentum spectra formed: 
in other words, how to describe a particle liber- 
ation process, which gradually destroys the local 
chemical and thermal equilibrium in expanding 
system. 



II. THERMALIZATION AND 
COLLECTIVE FLOWS 

Whether the interaction in the systems 
formed in A+A collisions are strong enough to 
result in a thermalization and collective effects 
such as hydrodynamic flows. A compatibility of 
the form of pion rapidity (y) spectra with that is 
predicted by hypothesis of the longitudinal hy- 
drodynamic flow was stressed first in the pio- 
neer Landau paper The interferometry sig- 
nature of the longitudinal flows, expressed in the 
specific rriT and y behavior of the long-radius: 
Rlong ~ J^l coshy at high p T [2], was con- 
formed by NA35/NA49 Collaborations (CERN) 
jj. The radial transverse flow in A+A collisions 
reveal itself in Nu Xu's plots showing a depen- 
dence of the transverse spectra slopes T e ff on 



hadron masses m (in the non relativistic approx- 
imation T e ff = Tf + where Tf M _ is the 
temperature and (v 2 ) is the mean squared trans- 
verse collective velocity at freeze-out). 

The most direct evidence of thermaliza- 
tion and transverse anisotropic (elliptic) flows at 
RHIC [J] is related to the behavior of the so- 
called V2 coefficients describing the anisotropy of 
transverse spectra in non-central A+A collisions. 
It is very naturally explained and basically de- 
scribed within hydrodynamic model for perfect 
fluid [5] wherein the initial geometrical anisotropy 
transforms into the collective velocity anisotropy 
of thermal hadronic gas due to different gradi- 
ents of the pressure in out- and in- directions as 
for the reaction plane at the initial stage. To 
reach a quantitative agreement with the experi- 
mental data one needs to use very small initial 
time, r ~ 0.5 fm/c, to start the hydro-evolution. 
If it starts at later times, neither the collective 
velocities, nor their (and the transverse spectra) 
anisotropy will be developed to a sufficient de- 
gree. 

These results, in fact, brought the two new 
ideas: first, that the quark-gluon plasma (QGP), 
at least at the temperatures not much higher T c , 
is the strongly coupled system - sQGP, and so 
behaviors as an almost perfect fluid, and, second, 
that thermalization happens at very early times 
of collisions. A necessity for the thermal pressure 
to be formed as early as possible appears also in 
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hydrodynamic description of the central collisions 
at RHIC. If one starts the hydrodynamic evolu- 
tion in "conventional time" t% =1 fm/c without 
transverse flow, the latter will not be developed 
enough to describe simultaneously the pion, kaon 
and proton spectra. The crucial problem is, how- 
ever, that even the most optimistic theoretical 
estimates give thermalization time 1 - 1.5 fm/c 
Q], typically, it is 2 - 3 fm/c. The discrepancy 
could be even more at LHC energies. 



III. PRE-THERMAL FLOWS, OR IS 
EARLY THERMALIZATION REALLY 
NEEDED? 

The way to solve the problem was proposed 
in year 2006, Ref. [7] and developed and exploited 
in Refs. [a, {J. It was shown that the initial 
transverse flows in thermal matter as well as their 
anisotropy, leading to asymmetry of the trans- 
verse momentum spectra, could be developed at 
the pre-thermal, either partonic/string or classi- 
cal field/ Glasma, stage with even more efficiency 
than in the case of very early perfect hydrody- 
namics. The illustration for the case of partonic 
free streaming is presented in Fig. [U 

The results of above mentioned papers show 

that: 

1) The radial and elliptic flows develop no 
matter whether a pressure already established. 
The general reason for them is an essential finite- 
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FIG. 1: The collective velocities developed in central 
(b=0) and non-central (b=6.3 fm) Au+Au collisions 
at the pre-thermal stage from proper time To = 0.3 
fm/c by supposed thermalization time n = 1 fm/c for 
scenarios of partonic free streaming in CGC model £<]/. 

ness of the system in transverse direction. Then 
the particle number or energy-momentum density 
flows directed outward the system cannot be com- 
pensated by the inward directed (from periphery 
to the center) flows. Just this creates non-zero net 
flows no matter how the collective velocity is de- 
fined: according to Ekkart or to Landau-Lifshitz. 

ii) The specific (linear) dependence of the 
transverse spectra slopes on the hadron masses 
in central A+A collisions and an anisotropy of 
the spectra in non-central ones happen only if, 
at least, partial thermalization happens. Then 
an almost isotropic particle emission in the local 
rest frames of the decoupling fluid elements gains 
Doppler blue shifts in the Lab system. These 
shifts depend on the collective velocities of fluid 
elements which are various in different directions 
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in non-central collisions - it leads to momentum 
spectra anisotropy. If no thermalization happen, 
similar as at free streaming, the final transverse 
momentum spectra will be practically the same 
as initial transversally isotropic momentum dis- 
tribution despite anisotropic collective flows in- 
deed develop in the non-equilibrated system. In 
that case also no specific dependence of the trans- 
verse spectra slopes on particle masses arises. 

Hi) So, the results, first published in 2006, 
show that whereas the assumption of (partial) 
thermalization in relativistic A + A collisions is 
really crucial to explain soft physics observables, 
the hypotheses of early thermalization at times 
less than 1 fm/c is not necessary. 



IV. PHENOMENOLOGICAL APPROACH 
TO THE PRE-THERMAL EVOLUTION 

Of course, the free streaming approximation 
for the processes at the very early pre-thermal 
stage is too rough and leads to a non local equi- 
librium structure of the energy-momentum tensor 
at the supposed "thermalization" time Tth- In 
the forthcoming publication [lo| we will present 
in details the phenomenological approach moti- 
vated by Boltzmann equations, which accounts 
for the energy and momentum conservation laws 
and contains the two parameters: supposed time 
of thermalization Tth and initial relaxation time 
T~rel( T o)- The result is: if some model or effective 



QCD theory bring us an energy-momentum ten- 
sor at the very initial time of collision tq, then 
we can estimate the flows and energy densities at 
expected time of thermalization Tth using equa- 
tions for the energy-momentum tensor of ideal 
fluid with (known) source terms which is associ- 
ated with the free evolving initial system when 
the interaction there is turned off: 



dA(x) = -Tf ree {x)d,V T ^ T {r) (1) 



with 



Vir) 



T f- T o 



Tf — T \ T rel(n)) 



tf - T 0, 

and T^ d correspond to T^ d of ideal fluid with 
renormalized energy density and pressure. 

% = - (1 " r(r))e,p - (1 - V(r))p). 



' hyd 



Such a method wil 



)e included in HydroKinetic 



Model (HKM) Q Q that accounts for non- 
equilibrium effects in similar way; the essence of 
the model is discussed in what follows. 



V. THE EVOLUTION IN THE 
EQUILIBRATED ZONE AND INITIAL 
CONDITIONS 

At the temperatures higher than the chemi- 
cal freeze-out temperature T c h the hydrodynamic 
evolution is related to locally equilibrated quark- 
gluon and hadron phases. The evolution is de- 
scribed by the conservation law for the energy- 
momentum tensor of perfect fluid: d v T^ u {x) = 
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and conservation laws for baryonic and strange 
net charge (Q) flows: d v Q v (x) = 0. In Ref. [11] 
the equation of state (EoS) at T > T c h is taken 
from Ref. [121] . that describes well the QCD lattice 
data at zero barionic chemical potentials and is 
matched with the chemically equilibrated multi- 
component hadron resonance gas at T = 175 
MeV. In this work we adjust that EoS to the 
parameters at the chemical freeze-out: [ib =29 
MeV, [i s =7 MeV, [i E =-1 MeV Q. 

For very central A+A collisions the initial 
transverse energy density profile for such an evo- 
lution is supposed to be the Gaussian as for the 
variable —r 2 /R 2 and with the maximal energy 
density e(r^,r = 0) = eo- We choose the ther- 
malization time to be Tth = 1 fm/c. The initial 
transverse rapidity profile, developed at the pre- 
thermal stage, is supposed to be linear in radius r: 
Ut = cej^. Note that the parameter a absorbs also 
a correction for underestimated resulting trans- 
verse flow since in this work we did not account 
for the viscosity effects neither at QGP stage 
nor at hadronic one. In the longitudinal direction 
the boost-invariance with the Bjorken flow is sup- 
posed. From the estimates presented in Ref. [9] 
R = 5.4 fm for Au+Au collisions. The only free 
parameters eo and a are defined from the fitting 
of the spectra and interferometry radii and turn 
out to be e = 17 GeV/fm 3 ; q=0.35. 



VI. THE EVOLUTION IN THE 
NON-EQUILIBRATED ZONE AND EOS 

At T < T c /j=165 MeV the system evolves 
as a chemically and thermally non-equilibrated 
hadronic matter. To calculate the EoS we 
approximate it below T c h by multi-component 
hadronic gas where to guarantee the correct parti- 
cle number ratios for all the quasi-stable particles 
a composition is changed only due to resonance 
decays into expanding fluid with possible recom- 
binations. We include 359 hadron states made of 
u, d, s quarks with masses up to 2.6 GeV. Thus, 
in addition to the equations of energy-momentum 
conservation, the equations accounting for the 
particle number conservation and resonance de- 
cays are added. If one neglects the thermal mo- 
tion of heavy resonances, the equations for parti- 
cle densities rii(x) take the form: 

d^{rn{x)u^{x)) = -Tim(x) + } j bjjTjnj(x) (3) 

3 

where bij = BijNij denote the average number 
of i-th particles coming from arbitrary decay of 
j-th resonance, Bij = Tij/Tjjot is branching ra- 
tio, Nij is a number of i-th particles produced in 
j — > i decay channel. The EoS in this chem- 
ically non-equilibrated system depends now on 
particle number densities rii of all the 359 par- 
ticle species. So the energy-momentum conser- 
vation equations and 359 equations ([3]) are solv- 
ing simultaneously with calculation of the EoS, 
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FIG. 2: The thermodynamic (e,p)- region occupied by 
the actual points (grey points) of equation of state 
(EoS) p(x) — p(e(x) , {ni(x)}) in chemically non- 
equilibrated multi- component hadronic gas during its 
evolution with IC as described in the body text. The 
different points are related to the various space-time 
points of J^-volume swept out by the expanding sys- 
tem. The lines corresponding to chemically equili- 
brated (red) and chemically frozen evolution (green) 
are marked for a comparison 

p(x) = p(e(x), {rii(x)}), at each point x. The 
thermodynamic region occupied by the EoS is 
presented in Fig. [5J Here the different points 
are related to the various space-time points of 4- 
volume swept out by the expanding system (with 
initial conditions as is described above). 

VII. SYSTEM'S DECOUPLING AND 
SPECTRA FORMATION IN 
HYDROKINETIC MODEL - HKM 

During the matter evolution, in fact, at T < 
T c h, hadrons continuously leave the system. Such 



a process is described in the HKM by means of 
the emission function S(x,p) which is expressed 
for pions through the term gain, G n (x,p), in 
Boltzmann equations and the escape probabilities 

CO 

Vtt{x,p) = exp(- / dsR 7r+ h(s±r+ - t),p)): 
S n (x,p) = G 7T (x,p)V 7T (x,p) [ill, lla]- For pion 
emission in relaxation time approximation G^ ~ 
f 7T R 7T+h -\ r GH^w where f n {x,p) is the pion phase- 
space Bose-Einstein distribution, R 7T+ h(x,p) is 
the total collision rate of the pion, carrying mo- 
mentum p, with all the hadrons h in the system 
in a vicinity of point x, the term Gn^n describes 
an income of the pions into phase-space point 
(x,p) due to the resonance decays. It is calculated 
according to the kinematics of decays with sim- 
plification that the spectral function of the res- 
onance H is 5{p 2 — (run) 2 )- The cross-sections 
in the hadronic gas, that determine via the col- 
lision rate R n +h the escape probabilities V(x,p) 
and emission function S(x,p), are calculated in 

n 

accordance with UrQMD method [16]. The spec- 
tra and correlation function are found from the 
emission function S in the standard way (see, e.g., 

a. 

VIII. THE RESULTS AND CONCLUSION 

The systems formed in A+A collisions go 
through different stages of evolution: from the 
initial one, which is far from equilibrium, to the 
thermal and chemically equilibrated phases of 
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sQGP and hadronic matter, then system again 
becomes non-equilibrated, neither chemically non 
thermally. At the late non-equilibrium stage 
the medium approximation for multi-component 
hadron gas is destroyed and hadrons are gradu- 
ally liberated. The correct description of both 
non-equilibrium stages: very initial one and the 
latest stage are very important since at the first 
the IC for the fireball evolution are generated and 
at the latest one the hadronic momentum spectra 
are formed. 

It is worthily note that at the first, 
pre-thermal stage the collective flows and its 
anisotropy are developed, as it is illustrated in 
Fig. [1] for the free streaming approximation. To 
go ahead one needs in the correct description 
of thermalization of the initially non-equilibrated 
system of either partons/strings or Glasma fields. 
The phenomenological approach dH) - ([21) [10] al- 
lows one to describe transformation from an arbi- 
trary initial energy-momentum tensor to the one 
corresponding to the perfect fluid in agreement 
with the conservation laws in the generalized re- 
laxation time approximation. 

The HydroKinetic Model allows one to de- 
scribe all the stages of the system evolution as 
well ELS EL formation of the particle momentum 
at the decoupling stage an agreement with the 
underlying transport equations. The basic hy- 
drokinetic code, proposed in is modified 

now to include decays of resonances into the 



expanding hadronic chemically non-equilibrated 
system. The non-equilibrium process of the spec- 
tra formation is evaluated in the first approxima- 
tion when the back reaction of non-equilibrium 
hadronic emission on the hydrodynamic evolution 
is ignored. In formalism of HKM l_l|] an account 



of the back reaction of the particle emission corre- 
sponds to an including of the viscosity effects into 
a description of hadronic hydro evolution. So here 
we estimate only the main characteristics of the 
dissipative effects in the most important region 
corresponding to the maximal pion emission. The 
results of our analysis show that for above men- 
tioned initial conditions for RHIC energies the 
maximal emission is achieved at the proper time 
r ~ 13 fm/c and transverse radii r ~ 8.4 fm. At 
the corresponding space-time vicinity, related to 
the last stage of the evolution where the spectra 



are formed, we estimate 



I/tt 



4.2 



J_ Vp 

Air ' s„ 



0.68 



Itt ' 



The ratio of the total (over all hadrons) shear 
viscosity to total entropy is ^ = 6.27^. The 
corresponding mean free pathes are = 3 fm, 
X p =1.2 fm and Knudsen numbers, defined as the 
ratio of the m.f.p. to the hydrodynamic length, 
are in this region Kn ff ~ 0.5, Kn p ~ 0.2. These 
estimates are important to understand the evolu- 
tion of multi-component hadron gas and applica- 
bility of (viscous) hydrodynamic approach to this 
stage. 

The crucial point for any model of space- 
time evolution is simultaneous description of 
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FIG. 3: Left panel: The transverse momentum spectra for negative pions, kaons and protons calculated in HKM 
for the cases of the vacuum cross-sections (solid lines), "transport" cross-sections (dash lines), with account 
for the resonance recombination processes (dotted line) and attractive for protons mean field (red line). The 
experimental data are taken from STAR [l9( and PHENIX [2l| Collaborations (RHIC BNL). Right panel: The 



HBT radii for the negative pion and negative+positive kaon pairs calculated in HKM for the same physical 
conditions as at the left panel. The experimental data: STAR [20( and PHENIX 



22 



23|. 



hadronic momentum spectra (in absolute units!), 
and measured space-time scales: the interferome- 
try, or HBT radii. The presented approach HKM 
leads to a good description of the experimental 
data as demonstrated in Fig. [3l Because of lack of 
the room, we demonstrate only R ou t/Rside ratios, 
- our results (not shown) for the transverse radii 
Rout and R s ide are also in good agreement with 



the data. As for the absolute values of the spec- 
tra, it is worthy to stress that parameters of the 
chemical freeze-out are chosen in accordance with 
STAR RHIC data, while for kaons and protons 
the PHENIX and STAR results are noticeably 
different. We analyze an influence of the main 
physical factors on the spectra and HBT radii. 
For this aim we consider interactions in multi- 
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hadronic gas with different cross-sections: vac- 
uum ones <Jij and with "transport" cross-sections, 
<Jij/2, and also we analyze the possibility when 
not only decays of resonances happen but also re- 
combination processes for non strange resonances 
takes place, we simulated the latter effect by 
utilization of the effective resonance widths in 
Eqs.flSJ), Ti — ► 0.75Fj, that corresponds to 50 per- 
cent increase of effective life-times of resonances 
due to the recombination processes. All these 
modifications shift significantly the time distri- 
bution of the hadronic emission but, as it fol- 
lows from Fig. [3j they do not influent essentially 
on pion and kaon spectra and the interferome- 
try radii. It illustrates and extend the general 
analytic results [17] that these quantities do not 
depend significantly on the freeze-out time at the 
isentropic evolution. As for the proton transverse 
spectrum, its best description requires, probably, 
the mean field contribution, which is attractive 
for protons [3] and so leads to reduction of their 
velocities in soft momentum region (central part 
of the fireball). The red dot line in Fig. [3J (left) 
corresponds to 14 percent reduction of the proton 
transverse rapidities in the region yx < 1. 

Summarizing we would like to emphasis that 
correct approaches to the non-equilibrated stages 
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